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ABSTRACT 

The presence of atomic gas mixed with molecular species in a "molecular" 
cloud may significantly affect its chemistry, the excitation of some species, and 
can serve as probe of the cloud's evolution. Cold neutral atomic hydrogen (HI) in 
molecular clouds is revealed by its self absorption of background galactic HI 21-cm 
emission. The properties of this gas can be investigated quantitatively through 
observation of HI Narrow Self- Absorption (HINSA). In this paper, we present a 
new technique for measuring atomic gas physical parameters from HINSA ob- 
servations that utilizes molecular tracers to guide the HINSA extraction. This 
technique offers a significant improvement in the precision with which HI column 
densities can be determined over previous methods, and it opens several new 
avenues of study of relevance to the field of star formation. 

Subject headings: astrochemistry — line: profiles — molecular processes — ISM: 
clouds — ISM: evolution 
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1. Introduction 

Star formation occurs in molecular clouds which are thought to have evolved from 
diffuse atomic hydrogen (HI) regions to form dense, cold, well-shielded regions composed 
primarily of molecular Hydrogen (H2). Our quantitative understanding of cloud evolution 
and specifically the conversion process (from HI to H2 gas) has been hindered by our inability 
to measure confidently the HI abundance in evolving clouds. In this paper we present a new 
technique for measuring HI column densities in dark clouds offering a significant improvement 
over previous methods. In the interests of brevity this paper includes the results from only a 
few clouds on which this technique has been applied in order to demonstrate the technique. 
The results of a much larger survey of observational data and analysis are to follow in a 
subsequent publication. 

The ability to determine accurately the HI component of molecular clouds could have 
a variety of benefits. Measurement of HI/H2 ratios in clouds, used in conjunction with as- 
trochemical models, allows us to determine the chemical ages of individual clouds or entire 
molecular complexes (e.g. Taurus, Perseus, etc.). This will greatly expand our understand- 
ing by constraining star formation models thus yielding insights into the collapse process 
and the interplay of magnetic fields, ambipolar diffusion, turblence, and various potential 
sources of cloud support. By studying age distributions in large scale regions we can learn 
about the processes that may trigger the collapse of large complexes. In contrast to previous 
methods, our technique allows us to determine the HI/H2 ratios for individual velocity com- 
ponents within a cloud thus yielding unique information about cloud kinematics. Further, 
the technique allows for the absolute measurement of quantities such as HI column density, 
in contrast to many previous studies which were limited to comparative measurements. 

HI is the dominant constituent of the diffuse ISM and its 21cm emission line is prevalent 
everywhere throughout the sky, especially near the galactic plane. Typical HI emission spec- 
tra are composed of numerous superimposed velocity components. The emission linewidths 
of the overall features are typically on the order of a few 10s of kms~^. They include velocity 
variations owing to galactic rotation as well as very significant peculiar velocities resulting 
from localized phenomena. It is rare to find a molecular region for which one can confidently 
claim that the HI emission observed along the fine of sight to the cloud is associated with the 
cloud, and is not due to background or foreground sources. Owing to such velocity crowding 
it is difficult or impossible to disentangle HI emission originating from a particular cloud 
from the background galactic HI emission. 

The situation for HI absorption is different. HI within a galactic cloud of any type 
can absorb the continuum emission from distant (galactic or extragalactic) radio sources. 
Because the HI optical depth varies inversely with the cloud temperature, the absorption by 
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galactic HI is stronger in cold, galactic HI clouds. The integrated optical depth of all the 
clouds along t h e line of si ght through the galac t ic dis k can exceed unity as demonstrated in 
Kolpak et al~l J2002h . and lCarwood fc Dickev I Jl989h . 



These cold, interstellar HI clouds may also be identified through their absorption of warm 
background HI emission originating within the galaxy. Because the emission being absorbed 
by cold HI clouds in this case is galactic HI emission, the resultant spectral absorption fea- 
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among others. It is important to emphasize that the term HISA refers 



to an observable spectral absorption feature rather than being a description of a specific 
physical process. 

Molecular spectral emission lines provide an independent view of that subset of inter- 
stellar clouds that are cold and composed primarily of molecular species. These "molecular 
clouds" are expected to maintain a residual abundance of atomic hydrogen if for no other 
reason than the cosmic ray disassociation of H2. This applies even when the chemical evolu - 
tion of the cloud has reached equilibrium (jSolomon k Werner Ill97ll : iGoldsmith k Li II2005I ). 
With the residual HI co-existing throughout the cloud with molecular species, the observed 
spatial and velocity (kinematic) structure of the molecular cloud will be similar whether 
observed in molecular emission or HI self-absorption lines. Thus, we expect to observe HI 
self-absorption features along the lines of sight of molecular clouds which share the spatial 
distribution and kinematics (non-thermal line width) of molecular emission lines. Such lo- 
calized association between i nolecular emission and HI se lf-absorption is observe d for many 
nearby clouds as reported by lLi k Goldsmith I (120031 ) and [Goldsmith k Li I (120051 ). The spe- 
cific case in which the HI absorption features observed in the direction of a molecular cloud 
share the spatial and kinematic structu re seen in the m o lecula r lines is called HI Narrow 
Self- Absorption (HINSA) as defined in iLi k Goldsmith I (120031 ). The term Narrow arises 
from the typically small nonthermal linewidths (on the order of 0.1 kms~^) of HINSA fea- 
tures, very similar to the similarly small nonthermal linewidths of molecular tracers along 
the same lines of sight. HINSA can be considered to be a subset of HISA, but it is a subset 
derived from an understanding of a specific, observable physical phenomenon in molecular 
clouds. 



While both are simply acronyms for spectral absorption features, HISA can be caused 
by a variety of different conditions and processes, but HI Narrow Self- Absorption (HINSA) is 
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that subset of HISA in which the atomic HI absorption correlates well with molecular emis- 
sion of certain tracers (most notably ^^CO) in sky position, central velocity, and nonthermal 
line width. Based on our current understanding of cold molecular clouds, the most satisfac- 
tory picture is that HINSA features are a result of HI gas located within these cold, dense, 
well-shielded regi ons. Some early examp le s of HINSA studies that pr edate the use of this 



term are those of IWilson fc Minn I (119771 ) , Ivan der Werf et al. I (119881 ) , and iJackson et al. 



(I2OO2I ). The technique which we describe in this paper pertains only to the extraction of HI 
data from HINSA features. 



The general picture which emerges, as found in lLi &: Goldsmith I (120031 ) and lGoldsmith Sz Li 



( l2005l ). is that the HI gas located within cold, quiescent cores of dark clouds produces HINSA 
absorption features. The well-defined center velocities and narrow line widths allow us to 
separate the HI gas associated with individual clouds from the galactic background. How- 
ever, the complexity of the background emission spectra that are frequently encountered 
makes extracting accurate data (especially in terms of obtaining the cold HI column den- 
sity) from the absorption features difficult. Several methods (discussed in fgj) have been 
used previously, but all are recognized to introduce significant uncertainties in the results. 
We here present a new technique that aims to improve the situation by using the properties 
of molecular emission to characterize the region producing the HINSA features, and then 
employes the HINSA spectral features to derive HI column densities. 

In ^ we present selected data and show the results of applying the new technique to 
them, and in ^ we contrast this with previous methods for analyzing HINSA data. In fJH 
and ^ we describe the technique and the combination of molecular data with the HINSA 
spectra. In ^ and fJT] we verify the validity of our technique using simulated data and also 
examine its limitations. 



2. Results Using Observational Data 

In this section we demonstrate the utility of our procedure by showing some results of 
its implementation based on observations of two molecular clouds. While a full discussion of 
the implications of these results is reserved for another publication, we hope that the results 
obtained encourage careful consideration of the technique described here. 

An observational survey, whose complete results and analysis are presented elsewhere, 
was performed using the Green Bank Telescope (GBT) and the Five College Radio Observa- 
tory (FRCAO) to obtain HI, OH, ^^CO, ^^CO, and C^^O maps of more than 30 dark clouds 
which exhibit HINSA features. Figure [T] shows sample HI spectra from two well known 
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clouds, Lynds 134 (L134) and L1757. HINSA absorption features are prominent in both 
sources. The red hues represent the recovered background spectra after HINSA has been 
removed. This representation is reassuring in that the smooth, natural looking background 
spectra indicate that the HINSA removal process has not produced marked distortions in the 
background. While L134 represents a simple case of a single emission component, in L1757 
our technique is able to discern two emission components closely spaced in velocity. Previous 
methods based solely on the HI line profile would have significantly overestimated the HI 
optical depth by assuming there was only a single emission component. By utilizing the 
associated molecular data as described below, the possibility of a single emission component 
is largely eliminated. 

The technique we use here assumes that the cold HI and molecular constituents are 
well-correlated spatially, so that we can assume that the line center velocity and nonther- 
mal line width for the cold atomic gas and the molecular gas in a particular direction are 
the same. This i s justi fied by previous studies of H INSA and CO isotopologues, including 



Li fc Goldsmith I (120031 ) and lGoldsmith fc Li I (120051 ). Our technique depends on use of these 



spectra to calculate the HI optical depths and column densities. In conjunction with ob- 
servations of molecular tracers this allows us to obtain HI column densities in dark clouds 
with improved accuracy. Further, we are able to obtain unprecedented detail by examining 
individual velocity components within each cloud, even when their velocity separations are 
small enough to make individual components difficult to distinguish in HI spectra. Addi- 
tional details are given in §1] 

Figures [21 [3l and Tables [1] and [2] show measured HI, H2, and total proton column 
densities, and the HI/H2 ratios for individual velocity components throughout numerous 
positions in clouds L134 and L1757. No such analysis has previously been carried out. The 
choice of these two sources indicates the importance and possible rewards of such studies 
in that the two clouds exhibit very different properties. While in L134 all the velocity 
components have similar HI/H2 abundances in the relatively narrow range between 10~^ 
and 10~^, in L1757 the situation is quite different, with ratios varying systematically as a 
function of velocity and extending over a much greater range. The most likely explanation 
for this phenomenon is that L134 is a mature cloud in which the chemical evolution of all of 
the velocity components has proceded to the point where HI abundances have approached 
their equilibrium values. In contrast, L1757 represents a much younger cloud in which each 
component is evolving at a different rate governed by its total density. This information 
yields valuable insight into the dynamics and chemistry of dark clouds, but has never before 
been available. A complete discussion of the results for a large sample of clouds will be 
presented in a forthcoming paper as a proper analysis requires more space than is available 
here. 
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L134 HINSA spectrum with solution 
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L1757 HINSA spectrum with solution 
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Fig. 1. — Sample HI spectra in the directions of the mlecular clouds L134 (upper) and 
L1757 (lower). We show the original observed spectra (black) and the recovered background 
spectra after HINSA removal (red), using the technique described in this paper. 
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L1 34 HINSA vs Total Proton Column Densities 
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Fig. 2. — Atomic to molecular abundance ratios in the molecular cloud L134. Each point 
represents a single velocity component at different position within the cloud. The 2.55-3.5, 
1.65-2.55, 1.55-1.65, 0.25-1.55, and -1.00-0.25 kms~^ velocity components are represented 
by black, red, blue, green, and cyan dots, respectively. These points are Nyquist sampled 
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L1757 Hl/Hj vs Total Proton Column Densities 
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Fig. 3. — Atomic and molecular hydrogen column densities and ratios in L1757, presented 
as in Figured The 4.3-6.7, 2.0-4.0, 1.0-2.5, -0.5-1.3 kms~^ components are represented by 
black, red, blue, and green points, respectively. Unlike L134, the various velocity components 
of L1757 show fairly well-defined but systematically different different relationships between 
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3. Limitations of Previous HI Self- Absorption Analysis Techniques 

Owing to the complexity of HI background emission along most lines of sight, it is clearly 
quite arbitrary to fit a straight line across the absorption feature. It is equally problematic 
to use OFF source observations since the HI background spectra change significantly over 
the angular size of the foreground cloud as sampled by the telescope beams that have been 
employed. Several methods for extracting HI column densities from HISA absorption have 
been tried in the past. These can be classified into two groups, each with its own set of lim- 
itations. The most commonly-used methods rely on fitting simple mathematical functions 
such as gaussians or pol ynomials over the observed spectra while masking the absorption 



features. As discussed in iLi fc Goldsmith I (120031 ) these techniques result in significant un 



certainty because the analysis is under-constrained. Arbitrary choices for the degree of the 
fitted polynomial and the masking ranges must be made, and these greatly affect the final 
results. 



Recently, iKavars et al. I (120051 ) implemented an automated algorithm for measuring 
HISA which relies on spotting the characteristic HISA features (the dip) in HI spectra. 
However, this t echni que (as pointed out by its authors) suffers from several limitations. The 



Kavars et al. I (120051 ) technique, as well as function fitting methods, are unable to deter- 
mine the gas temperature and the optical depth without assuming a value for one of these 
important quantities. These methods are thus limited to comparative measurement of one 
quantity, such as HI column density, over large regions recognizing that absolute values 
cannot be obtained. 

Owing to the complexity of the galactic HI background, and lacking additional in- 
formation, it is often impossible to tell which dips in the spectra are genuine absorption 
features and which are produced simply by the superposition of two emission features. Both 
phenomena may be spatially extended with continuous velocity distributions. Molecular ob- 
servations can be used to distinguish between these two scenarios. In molecular clouds there 
are often multiple components along the same line of sight with closely spaced central veloc- 
ities. The lower mass of atomic hydrogen compared to molecular tracers (and consequently 
greater thermal line width) means that it is often possible for multiple components which 
are distinguishable in molecular spectra to merge into a single HI absorption feature. Using 
the molecular g 8bS clS 8b template for the HINSA analysis has the potential to greatly en- 
hance the accuracy of the HI/H2 comparison since the ratio may be different in each velocity 
component. 

Previous analysis techniques have been applied to both HISA and HINSA sources with- 
out distinguishing between the two. Our technique is limited to HINSA sources for which 
(by definition) there is molecular data available. Limiting our scope in this manner allows 
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us to utilize observational data from various molecular tracers to assist in extracting mean- 
ingful results from HINSA features. As will be shown in what follows, to the extent that our 
assumption that HI is mixed with the molecular gas is valid, we are able to measure both 
the temperature and optical depth of the HI gas without any ambiguities regarding fitting 
functions. In this way we may treat each velocity component in a cloud separately. 

4. The Technique 
4.1. Analytic Representation 

An idealized HINSA spectrum can be represented by the radiation transfer equation as 

TAiv) = iniv) + T,)e-<^^ +Th{1- e'^^^)) , (1) 

where v corresponds to velocity, Ta{v) is the observed spectrum, T^iv) is the background 
emission due HI emission clouds, represents the continuum emission along the line of 
sight at 21 cm including the 2.7K microwave background ( §4 . 2 p . Th is the gas temperature 
of the foreground absorbing cloud producing HINSA, and t{v) is the optical depth of the 

absorbing HI gas at 21cm which we can describe as t{v) = tqc '^'^h ^ where tq is the peak 
optical depth, vh is the central velocity of the HINSA component, and an is the HINSA 
linewidth. In this particular representation we have assumed that there is no intervening 
HI emission originating between us and the absorbing HINSA cloud; this topic is further 
discussed in §4.3[ When multiple absorption components overlap in velocity within the same 
line of sight, the transfer equation must be applied consecutively for each component from 
farthest to nearest. 

HINSA absorption features are typically gaussian-like dips with much smaller ampli- 
tudes and linewidths than the characteristics of the emission spectra (Tfe(f)) on which they 
reside. These features can be so small in comparison to the background emission that they 
can be difficult to distinguish by eye. However, by utilizing a fundamental property of gaus- 
sians we can obtain a much better representation of the data in which the HINSA becomes 
more apparent. 

The derivatives in velocity or frequency of a gaussian function can be written as: 

fiv) = he^-"'"-'^ ^ = -^f{v) ^ = -j,f{v) + ^^f{v) , (2) 
where it is apparent that higher order derivatives show a stronger dependence on 
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linewidth. For conceptual purposes we can temporarily approximate the HINSA absorption 
and background emission features as gaussians. Because the linewidths for HINSA features 
(typically 1 kms~^) are much narrower than the linewidths for the background emission 
components, in the first and second derivative representations, even a low-amplitude absorp- 
tion feature will become readily apparent in comparison to the bright background emission. 
This is illustrated in Figure HI In the second derivative representation, HINSA becomes the 
dominant feature in the spectrum thus making it easier to isolate and quantify. 

One could in principle use higher derivative representations to search for the best fit as 
the HINSA becomes progressively more prominent, however with real data the noise also 
rises with higher derivatives. The second derivative representation has been found to be 
the optimum choice for our GET and Arecibo observations. It would be possible to take 
the second derivative representation of any HINSA spectrum and attempt to fit a simple 
function to it to extract the absorption. However, HINSA spectra are rarely so simple as to 
lend themselves to this kind of fitting and we would have very limited information on the 
temperature of the HI gas. Furthermore, the large thermal linewidths of HI spectra tend 
to blend all the velocity components along a single line of sight. Unless we are convinced 
that all the velocity components in a cloud can be characterized by the same HI/H2 ratio, 
we must find a way to extract the HINSA using additional constraints based on physical 
parameters. 

The primary benefit of limiting our technique to HINSA features is that we can use 
data from molecular tracers in order to constrain the analysis. If we can determine the 
values of all the other parameters (linewidth, center velocity, and temperature) for each 
velocity component then we can ideally limit ourselves to searching for that value of tq which 
produces the emission spectrum whose second derivative best removes the HINSA feature. 
Minimizing the integrated squared sum of the second derivative background spectrum is 
sufficient to extract the HINSA feature (see §4.51) . 

4.2. Continuum Emission at HI wavelengths 

There are various sources of continuum emission (Tc) in equation [1] throughout the 
galaxy that contribute to the sky brightness at 21cm. The emission may originate from a 
variety of sources such as synchroton emission within the galaxy and continuum emission 
from ionized gas such as that in the Galactic Halo. Continuum emission can also originate 
from individual discrete galactic and extra-galactic sources, such as pulsars, HII regions, 
AGNs, accretion disks, and quasars. These sources will typically produce emission whose 
intensity is constant across the observed HI band. Observing techniques such as frequency 
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An idealized HINSA spectram and its derivatives 
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Fig. 4. — (upper) An idealized HINSA spectrum with typical observed properties along with 
its first and second velocity derivatives. The scales have been normalized for demonstration 
purposes, (lower) Second derivative representation of the observed spectrum toward L134 
shown is Figure [H It is apparent that a small HINSA dip becomes the dominant feature in 
the second derivative representation 
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switching as well as software baseline removal algorithms remove the continuum emission 
and the 2.7K background along a particular line of sight from the reduced spectrum. The 
continuum emission component should ideally be included when calculating HINSA since the 
temperature of the cold gas clouds ( lOK) is comparable to the antenna temperature pro- 
duced by the continuum emission. Several galactic sur veys have been perf or med which have 



measu red t he value of the 21cm c ontinuum including iReich et al. I (119971 ). lUyaniker et al. 



(I1999I ). and iReich fc Reich I (119861 ). among others. These data can be utilized to determine 
the expected continuum emission along any HINSA source. Typical values are on the order 
of a few K. According to equation [1], the value of Tc, while significant, becomes critical only 
when the HINSA optical depths are high and gas the temperature is low. 



4.3. Foreground HI Emission 

Usually there will be HI emission between the observer and the HINSA cloud due to 
intervening HI gas. To assess the importance of the effect this invervening gas might have on 
our HINSA measurements we must expand equation [1] to include all three gas components, 
that of the HI emission behind the HINSA cloud, absorption due to the HINSA cloud itself, 
and emission between us and the HINSA cloud. 



Tj{v) = T.e-^"'^^^ + T,, b(l - e-^"'^^^) (3) 

represents the emission component due to the warm HI gas behind the HINSA cloud whose 
excitation temperature is Tex b and optical depth is Tb{v). At this stage, only the background 
continuum emission is attenuated. 

The HINSA cloud then attenuates the background emission, while also adding its own 
contribution as in 

Tu{v) = Ti{v)e~^"^^^ + Tex h{1 - e'^^^")) , (4) 

where th and Texn represent the optical depth and excitation temperature of the cold HI 
gas within a molecular cloud. Similarly, any foreground gas will attenuate the spectrum 
further as well as adding its own emission to produce the final observed spectrum: 

Tni{v) = Ta{v) = Tn{v)e--f^-^ + T^, ^(1 - e~^f^^^) , (5) 

where Tf and T^x / represent the optical depth and excitation temperature of the foreground 
HI gas between us and HINSA cloud. 

We now have two new terms, Tf and T^x /, describing emission due to HI gas between 
us and the HINSA cloud. There is no way to operationally disentangle which portion of 
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the observed emission spectrum originates from in front or behind the molecular cloud, but 
clearly any foreground emission we observe will attenuate the HINSA features and we must 
somehow correct for this. To get a sense of the effect and scale of the attenuation, we can 
assume that both the background and foreground emission gas have the same excitation 
temperature Tex tot = Tex b = Tex /• We can further assume the HI gas is distributed 
uniformly throughout the galaxy along the line of sight so that a linear relationship can be 
estabhshed between and Tf such that = pr^ot and tj = (1 — p)Ttot where r^f is the 
total optical depth of all the emitting HI gas behind and in front of the HINSA cloud and p 
represents the fractional distribution between the two which we can estimate using 



disk 



Here, Ddoud is the distance to the HINSA cloud, b is the galactic longitude along the given 
line of sight, and Rdisk represents the distance from Earth to the edge of the galactic disk 
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a more 



we are using geometry of the local galactic disk, and the HINSA cloud's position within the 
disk to estimate how much of the HI emission visible along a particular line of sight is in 
fact located in front of a HINSA source attenuating the HINSA feature. 

This emission, if located at the same velocity as the HINSA feature, will have the effect 
of washing out the absorption feature and thus the apparent HINSA optical depth will be 
lower than the intrinsic value. In fact, it can be s hown from the above formu lation that for 
optically thin HINSA, p{l,b) oc r^^ as showin in Ivan der Werf et al. I (Il988l ). where r now 
represents the observed apparent optical depth of any HINSA feature. Thus to obtain the 
intrinsic HINSA optical depth of a cloud we must divide the apparent optical depth by p{l,b). 
Due to the uncertainties in the HI distribution it is impractical to directly model the effects 
of foreground clouds when trying to extract HINSA profiles. Instead HINSA optical depths 
are computed on the assumption that there is no intervening emission and are subsequently 
adjusted using the factor p{l,b) which itself depends only on the cloud's position and is 
calculated independently of the observed emission spectrum. 

The observed spectra thus give us an apparent optical depth which then must be ad- 
justed by p to yield an estimate of the intrinsic optical depth. For most nearby clouds this 
correction factor is in fact small in comparison to other sources or error. It must be true that 
p{l,b) is not much less than unity if a significant HINSA feature is to be observed otherwise 
the foreground emission would obscure the HINSA absorption. For the two sample clouds 
in this paper, L134 and L1757 the correction factor has values of 0.5 and 0.7 accordingly. 
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4.4. Obtaining Constraints Using Molecular Parameters 

It is possible to obtain data from molecular tracers which can provide us with reliable 
estimates of the HINSA feature proper ties for each veloc i ty com ponent by utilizing certain 



assumptions justified by previous work. iLi &: Goldsmith I (|2003| ) showed a very strong corre- 



lation between the line center velocity of each HINSA component and those of its molecular 
counterparts using the OH, ^^CO, and C^'^O molecules as tracers. Thus, once each molecular 
component has been identified and separated using traditional fitting we can assume that 
each velocity component visible in molecular emission has an accompanying HINSA compo- 
nent whose optical depth may yet turn out to be 0. We can use the central velocity from 
any of the above molecules as an estimate of the number of HINSA components present and 
their center velocities (vh)- 

Spectral linewidths are composed of both thermal and nonthermal(turbulent) compo- 
nents. The two components add according to aobs = a/ crff^ + cr^^, where aobs, crth, and ant are 
the observed linewidth and the thermal and nonthermal components respectively. The ther- 
mal component arises from the random motions of individual particles within the gas and 
will be different for varying atomic or molecular tracers due to the different particle masses. 
The nonthermal component results from turbulence and bulk motions within the cloud. If 
two species of gas were uniformly mixed throughout a cloud one would expect their spectra 



to share similar non-thermal linewidths. iLi fc Goldsmith I (120031 ) showed a good correlation 



between the nonthermal linewidths of HINSA and the accompanying ^^CO emission. We can 
use ^^CO to estimate our expected HI linewidths, but first we must disentangle the thermal 
components using information about the gas temperature within the cloud. 

Under the assumption that the cold HI is in thermal equilibrium with the surrounding 
molecular gas we can use the kinetic temperature of a molecular tracer to estimate the HI gas 
temperature (Th)- The J = 1 — transition of ^^CO has been used in the past to estimate 
gas temperatures because it usually appears as optically thick in molecular clouds. Equation 
[T] is valid for molecular emission as well as HI absorption as it describes generic radiative 
transfer. In the case of ^^CO there is generally no background or continuum emission and 
so the equation reduces to 

TA{v)=Ti2co{l-e-^''coi-)) , (7) 

where Ta{v), Ti2co, and ti2co{v) describe the observed ^^CO antenna temperature, the 
^^CO gas excitation temperature (which in high density regions can be taken to equal the 
kinetic temperature), and the optical depth. If the optical depth is taken as ^ 1 then the 
equation reduces simply to T4 = Ti2(jq. Thus ^^CO provides us with estimates of the HI 
gas temperature allowing us to estimate the HINSA optical depths, column densities, and 
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thermal and non-thermal linewidths. Some derived ^^CO temperatures for L134 and L1757 
are given in Table [2l Since ^^CO emission has been shown to exist wherever there is ^^CO 
emission and HINSA we can measure the ^^CO temperature at each position and use it 
to estimate the ^^CO and HINSA temperatures at each position in the cloud. Using these 
constraints we are left with only one free parameter, the peak optical depth of the cold HI 
gas. 

Due to the uncertainties in the assumptions made the constraint parameters cannot be 
taken as exact values of the HI gas properties. Instead they can be used as close guess values 
to help us constrain the fit. By using ^'^CO and ^^CO data we are also able to estimate H2 
column densities using fractio nal abundances of carb on monoxide that are well-determined 



in many nearby dense clouds (IStahler fc Palla II2004I ). 



4.5. The Fitting Procedure 

As explained in section H] any narrow HINSA feature will be the dominant feature in the 
second derivative representation of the observed spectrum. As such, once guess parameters 
for the center velocity, line-width and temperature of the HI gas have been obtained from 
molecular tracers the remaining task is to find the value of tq for each velocity component 
which minimizes the resulting area under the curve of the second derivative of the recovered 
background emission spectrum after HINSA is removed in accordance with the discussion 
in §4.1[ This is done by first using the trial values of the HINSA linewidth, center velocity, 
temperature, and optical depth and using those on a particular observed spectrum (Ta{v)) to 
obtain the recovered background spectrum without HINSA (Ti,{v)) as described in Equation 
dJ According to §4.11 the HINSA feature will be revealed as a large feature in the second 
derivative representation of any observed spectrum. To find the best trial Tb{v) spectrum the 
algorithm must search for one whose integraded intensity in the second derivative represen- 
tation is minimized. The second derivative of the recovered spectrum {Th{v)) is calculated 
numerically through traditional methods. Since the second derivative representation will 
have negative values as well as positive ones, it is squared prior to integration. We thus 
minimize the value of / as described by the function 




which is simply the squared integrated intensity of the second derivative of the background 
spectrum (T^) that is recovered after trial values of the HINSA optical depths and other 
parameters are removed from the observed spectrum. Finally, determining the smallest 
value of / will give us the best-fit values of the HINSA parameters. It is possible to use a 
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more sophisticated minimizing function, but the HINSA feature is so dominant in the second 
derivative representation that this has proven not to be necessary. 

Previous studies have shown that the correlation between HI and ^^CO center velocities 



is strong (ILi fc Goldsmith II2003I ). Hence, the vh parameter is held fixed for each velocity 
component. The Tp estimates obtained from ^^CO do not carry as much confidence (as 
discussed in the next section), however the temperature is still kept as a fixed parameter 
during our fit. For optically thin spectral lines such as HINSA the temperature cannot 
sufficiently constrain the shape of the absorption feature to be allowed as a free parameter 
since a wide range of combinations of temperatures and optical depths would be able to fit 
the data. Thus the gas temperature must be specified in any fits. 

The linewidths obtained from ^^CO observations can only be applied to the HI gas under 
the assumption that the two species are coextensive throughout the cloud. This assumption, 
while strongly suggested by the previously-referenced studies, is clearly not exact. Hence 
the molecular linewidths for each velocity component are used as guess values for the HI fit. 
The HI linewidths are allowed to vary in value during the fit near the initial guess values. 

In practice, tq is left as a free parameter, and an as a constrained parameter for each 
component. For every trial value of tq a T;, spectrum is constructed and its second derivative 
is checked for the fitting parameter as specified by equation [HI All velocity components 
along the same line of sight must be fitted simultaneously to account for obscuration due to 
overlapping. There are a variety of fitting algorithms that can be used for the final fitting. 
In the present case a variation of the Levenberg-Marquardt method was implemented. In 
the interest of clarity, we give in appendix [A] a step-by-step description of the entire process 
including the determination of molecular parameters and HINSA extraction. 



5. Coming to Terms with Variable Beam Sizes and Unresolved Sources 

Our technique combines observations from different telescopes with different frequencies 
and beam sizes. The HI data presented in this paper was obtained using the 100 m Greenbank 
Telescope (GET) having a 9' beam size, while the ^^CO and ^^CO data were obtained using 
the Five College Radio Observatory (FCRAO) 14 m telescope having a 45" beam size. 
Combining these data sets is non-trivial. Convolving the CO data directly to a 9' beam has 
the effect of smoothing out a great deal of structure, kinematically as well as spatially, which 
is neither necessary nor desirable. 

By fitting the molecular data prior to convolution we preserve information regarding 
individual velocity components which may otherwise disappear or be distorted in an HI beam. 
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It is not necessary for all molecular components to have accompanying HINSA features. 
Maintaining this information is crucial to studying the comparative composition of different 
velocity components within a cloud as well as properly determining their HI/H2 ratios. High- 
resolution molecular data should only be convolved to a beam size which results in sufficient 
signal-to-noise to produce good fits to the individual emission components. Following this 
argument and a series of numerical experiments, our FCRAO spectra were convolved to a 
2 arcminute beam size. Parameters including center velocity, linewidth, temperature, and 
optical depth for each molecular emission velocity component are derived using conventional 
fitting techniques. Corresponding guess parameters for the HI gas are then constructed by 
convolving the molecular fit parameters to the HI beam. This approach allows us to measure 
the HINSA column densities for each individual molecular component, as well as to utilize 
the high spatial resolution information on the molecular constituents of the cloud. 

Many s ources are unresolved w hen o bserved with available s ingle-dish HI beams. Previ- 



ous studies ( iLi fc Goldsmith I (120031 ) . and lGoldsmith fc Li I (120051 )) have shown strong spatial 



and linewidth correlations between HINSA and ^^CO. With this information we may make 
the assumption that HINSA and ^^CO are well mixed throughout a cloud. This affords us 
the opportunity to make estimates of the HI beam filling factors by assuming that HINSA is 
present wherever significant ^^CO emission is observed. This assumption, while not precise, 
is useful. In the case of resolved sources the beam filling factor plays a minimal role. In 
unresolved sources, or the edges of larger clouds, the beam filling factor as implemented by 
the present technique, does not alter the calculated HI/H2 ratio, but only alters the total 
column densities of both species. 



6. Example Solutions Using Simulated Data 

A controlled way to assess the validity of our technique is to apply it to synthetic, yet 
reahstic data and compare the derived results with the known input values. Simulated data 
does not allow us to test the assumptions on the relations between molecular and atomic gas 
that are used by this technique. It is possible however to assume those relations to be correct 
and to test the validity of the HINSA fitting technique by itself. Further it is possible to 
observe the effects on the HINSA fit of having incorrectly determined the fitting parameters 
based on molecular data. 
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6.1. A Simple Case 

The simplest case consists of a single emission component obscured by a single HINSA 
absorption feature as in Figure [51 As expected, HINSA is the dominant feature in the second 
derivative representation. For the simulated run, the center velocity and temperature of the 
cold HI gas are provided as precisely known quantities. The linewidth, however, is given 
only as a guess parameter. The task then is to search for those values of the linewidth and 
the peak optical depth which minimize the squared sum of the second derivative background 
spectrum after HINSA is removed. Figure [6] shows the original observed spectrum as well as 
the recovered background spectrum after fitting. The resulting computed optical depth of 
0.198 compares favorably to the original inputed value of 0.2. It is apparent from the second 
derivative representation of the recovered spectrum that the HINSA which was prominent 
in the original spectrum (Figure [5]) has been effectively removed. 



6.2. Effects of Faulty Parameter Determination from Molecular Data 

It is clear from the previous simulation that in an idealized case where the properties 
(center velocity, linewidth, and temperature) of the HI gas are well known, the HINSA optical 
depths can readily be determined to a high degree of accuracy. However in practice, the 
properties of the HI gas are not known but are instead surmised with some uncertainty from 
observations of molecular tracers. Hence it is necessary to examine the impact of incorrectly 
deriving the HI properties from molecular data. One concern is the gas temperature, which 
we obtain from ^^CO. It is reasonable to assume that the HI excitation temperature is equal 
to the kinetic temperature of the gas. However, temperature gradients may exist resulting 
in the average temperature appropriate for analyzing the HINSA being different from that 
determined from the optically thick common isotopologue of carbon monoxide. Since HINSA 
features are typically optically thin, an incorrect estimate of the HI gas temperature would 
not strongly alter the shape of the fitted HINSA, but would only affect the derived optical 
depths and column densities. Thus, incorrect temperature estimates would be difficult to 
detect through fitting alone. 

Incorrect estimates of the HI linewidths based on ^"^CO data are plausible since the 
relation between them is only approximate. Here we investigate the effect of such an error 
by forcing our fitting routine to find the best fit HINSA optical dep t h whi le being constrained 



to use an incorrect value for the HI linewidth. iLi fc Goldsmith I (120031 ) showed a variance 
of roughly 20% in the observed relation between HINSA and ^'^CO nonthermal linewidths. 
Figure [7] shows the resulting fit when forced to utilize a HINSA linewidth 20% larger than 
the actual value of the same simulated data as in the previous example. It is no longer 
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Fig. 5. — An idealized HINSA spectrum (red) from simulated data. The blue line represents 
the background HI emission spectrum, while the green line represents the second velocity 
derivative of the observed HINSA spectrum. The spectra include typical radiometric noise. 
The HINSA feature has a peak optical depth of 0.2 and gas temperature of lOK. 



l_ 

O 
i_ 



O 
C 

C 
< 



100 



50 



Observed Spectrum with HINSA 
Recovered Background Spectrum after HINSA removal 
Second Derivative of Recovered Bacliground Spectrum 




Inputed HINSA Optical Depth: 0.2 
Derived HINSA Optical Depth: 0.1 9S 



-50 



-20 



■10 



10 



LSR Velocity (km/sec) 



20 



Fig. 6. — The original observed spectrum (blue) along with the derived background spec- 
trum after HINSA removal (red) using the same simulated data as in Figure [51 The recovered 
background spectrum, and the calculated optical depth of 0.198 are nearly identical to the in- 
put parameters. The second derivative representation of the recovered background spectrum 
(green) shows no signs of the HINSA feature. 
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possible for the algorithm to produce a smooth recovered background spectrum, regardless 
of the optical depth assumed. This is a key indicator of an error in the input information. 
Figure E] shows that the new best fit HINSA optical depth is 0.245 rather than 0.20. Similarly 
inaccurate solutions result when a faulty center velocity is used. 

6.3. Complex Examples 

It is often the case that two emission components may overlap at the velocities at 
which HINSA may be expected, thus obscuring the absorption feature and increasing the 
uncertainty in the analysis. In such cases, previous techniques have had great difficulty in 
determining whether any HINSA is present at all. Our technique attempts to recover the 
smoothest, most gaussian-like background spectrum while being constrained by the HI gas 
parameters derived from molecular observations. It is unlikely for emission components to 
appear over a HINSA features in such a way as to completely obscure the absorption. 

As an example. Figure [9] shows that our technique is able to accurately extract HINSA 
features which on visual inspection might be ambiguous. In this simulation, the HINSA fea- 
ture with an optical depth of 0.20 was placed in between two identical emission components. 
The recovered HINSA optical depth had a value of 0.208 which is well within the error limits 
of other uncertainties as discussed in the section. Figure [TD] shows further examples with 
more complex geometry in which our technique is still able to pick out HINSA features with 
good precision. 

6.4. Distorted Emission Components 

The technique we present here is based on the assumption that the background emission 
is composed primarily of a superposition of gaussian-like components. Significant deviations 
from this assumption can lead to incorrect results. Figure [11] shows the case in which the 
background spectrum is complex, due to the superposition of overlapping components. The 
presence of a background component of width only modestly greater than that of the HINSA 
feature, in combination with the broader emission feature, results in a noticeable error in the 
derived HINSA optical depth. On casual inspection, there may be little indication that the 
solution is problematic (as shown in the lower panel of Figure [T2l) . When dealing with real 
data, where the actual background spectrum is not known beforehand, there is no simple way 
to determine when such distortions may occur, or what their effect might be on our derived 
HINSA optical depths. Though instances such as this are likely to be rare, it is worthwhile 
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Fig. 7. — Using a HINSA linewidth 20% greater than the correct value when fitting the same 
simulated HINSA spectrum as in Figure [5] (blue) produces a distorted recovered background 
spectrum (red). The second derivative representation has a clear signature at the velocity of 
the HINSA absorption feature observed, and both its apearance and that of the the recovered 
background spectrum are unhkely to be representative of a realistic system. 
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Fig. 8. — The squared sum residuals of the second derivative of the recovered background 
spectra for different trial values of the optical depth using the same data, and incorrect 
linewidth as in Figure [71 This function I represents the quantity to be minimized in our 
fitting procedure described in §4.5[ The minimum value of I is considered to be the best fit, 
and in this case corresponds to an optical depth of 0.245, compared to the correct optical 
depth of 0.20. Realistic simulated statistical noise is included in this plot, and while the range 
in optical depth corresponding to a specified increase in the residuals may be significant, there 
are no localized minima to confuse the results. 
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Fig. 9. — Two relatively broad identical emission components combined with a HINSA 
feature located midway between them yield the simulated observed spectrum (blue). The 
recovered background spectrum (red) and accompanying second derivative representation 
(green) illustrate the results after HINSA extraction. The derived optical depth is 0.208, 
which compares very well with the input optical depth of 0.20. 
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Fig. 10. — Upper: Two unequal emission components together with a 0.20 peak optical 
depth HINSA feature result in the simulated observed spectrum (blue). The recovered 
HINSA optical depth is 0.215, and the HINSA extraction results in the recovered background 
spectrum (red), which is nearly identical to the input background spectrum. Lower: A similar 
case where the HINSA feature has been placed to coincide with a narrow emission component 
whose shape is slightly distorted by the presence of another emission component. The derived 
HINSA optical depth of 0.232 is still close to the value 0.20 that was input to the calculation. 
The second derivative representation of the recovered background spectra (green) shows that 
there is still significant structure left after HINSA extraction. This structure is the result of 
the interaction between two emission components. 
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to be especially wary of HINSA measurements when the absorption feature seems to lie in 
the trough between two emission components. 

6.5. Confidence in Derived HINSA Optical Depths from Simulated Data 

It may be tempting to take Figure [S] and use standard methods for determining the 
uncertainties in least-squares linear fitting problems to determine the uncertainties in our 
estimates of r based on the width of the trough of the fit residuals. But our problem does 
not strictly allow such analyses. There are two components that contribute uncertainty 
to our estimates of r. The first is radiometric noise. The simulated data in this section 
all include realistic simulated noise for observations made with a large single-dish telescope 
such as Arecibo or the GBT. This source of uncertainty is in fact linear in nature and could 
be treated using traditional least-squares linear fitting analyses. However, as HI spectra are 
typically bright and easily obtainable with only a few seconds of integration time this source 
of uncertainty is insignificant compared to the other uncertainties. If we were to look directly 
at Figure [H then the uncertainties due to radiometric noise would appear as small, narrow 
localized minima. However, in the case of this and all our simulations this contribution is 
smaller than even the thickness of the line in the plot. 

The dominant source of uncertainty which produces the trough is non-linear and some- 
what systematic. The fundamental principle on which our technique depends, as discussed 
in §4.11 requires that our observed spectra resemble the superposition of multiple (positive 
and negative) gaussians. According to the radiative transfer equation however this is not 
strictly the case. When multiple emission and absorption components are superimposed they 
attenuate one another in a non-linear fashion. Hence our observed spectra will never strictly 
represent gaussians, but will only resemble them. This introduces a somewhat systematic, 
" nonlinearity" error which is much greater than that introduced by radiometric noise. 

As an example, if we were to re-run the simulation from Figure [H] numerous times with 
all the same parameters for emission and absorption components but different radiometric 
noise (of the same average amplitude), the measured optical depth would always be very 
close to 0.245, the only variation in the results would be due to radiometric noise which 
in these realistic cases is insignificant. This represents a 22.5% error (not an uncertainty) 
relative to the original inputed value of 0.20 which is produced not by random local minima 
which would confuse the results, but by a systematic shifting of the entire residuals trough. 
This systematic error is produced when several emission and absorption components are su- 
perimposed in such a way that the final spectra no longer resemble simple superpositions of 
gaussians. Throughout this section we have shown successively more complex superpositions 
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Fig. 11. — Upper: For comparison with Figure [IHl a third emission component with bright- 
ness temperature of 30 K and center velocity of kms~^ has been introduced to distort the 
shape of the emission spectrum at the location of the simulated HINSA feature. In this 
instance the derived HINSA optical depth of 0.145 differs considerably from the inputed 
value of 0.2. Lower: An expanded view of the region of interest showing the recovered and 
inputed background emission spectra. 
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of this nature in an attempt to test the hmits of this systematic error. We have found that 
for the vast majority of observations hkely to be encountered, this error will be only a few 
percent. In an attempt to find the limits, we constructed the most complex superpositions 
that would likely be observed in real data and found that in the very worst cases the system- 
atic error reached only about 50%. Figures [H [TOl and [11] represent some of those complex 
cases. For the majority of observations, this nonlinearity error will be minor in comparison 
to other sources of error due to uncertainties in physical conditions in the cloud. In the most 
extreme cases, the nonlinearity error will be on a par with, but not dominant compared to 
these other sources of uncertainty, which are discussed in the following section. 

These simulations cannot directly test the validity of our assumptions regarding the 
relations between atomic and molecular gas properties as those are inherently non-statistical 
uncertainties whose properties are currently not well known. However, the simulations do 
show that the solutions are modestly sensitive to any such errors in that errors in our ^^CO- 
based estimes of either the HI linewidths or center velocities neccessarily produce visible 
artifacts in the recovered background spectra after HINSA extraction, as shown in Figure [71 
This allows us to identify faulty fits, where we can then go back to the molecular data and 
construct better templates for the HINSA extraction. Hence it is important to understand 
all the possible sources of uncertainty in our estimates of the HINSA properties based on 
molecular templates as discussed in the following section. 



7. Limitations: Sources of Ambiguity and Uncertainty 

While it has been shown in the previous section that our technique is able to produce 
good fits to the HI data when the HI gas properties are correctly derived, there are several 
factors which limit the accuracy of our technique. These sources of error are non-statistical 
in nature in that they arise from our reliance on the association between HI and molecu- 
lar tracers. As such, they make the placement of simple error bars on derived HI optical 
depths impossible. While it is certainly possible to place uncertainty estimates numerically 
through standard least-squares fitting methods, and compound those with estimates of the 
uncertainties due to measurements of the molecular parameters and their relation to the 
HI gas, any such numbers would be meaningless and misleading. The best way to estimate 
the uncertainties in our results is to examine the spread in the results in Figures [2] and [3] 
and assume that in an ideal case each component would be represented by a single line (of 
certain shape) which would be governed by the physical processes within the cloud such as 
the HI to H2 conversion process. This would give us an upper limit on our uncertainty in 
deriving the HI/H2 ratio. However, the present models of cloud evolution do not yet offer 
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us sufficient confidence to make sudi estimates witfi certainty. More detailed observations 
and modeling of molecular cloud processes are necessary before such estimates can be made 
with confidence. It is our hope that the technique described in this paper will help in this 
regard. This development is likely to be an iterative process. It is important to note that 
the sources of uncertainty in previous techniques, while inherently somewhat different, are 
not smaller nor better determined than those in the method presented here. It is, however, 
of critical importance to quantify and to understand each of these sources of uncertainty in 
order to understand the level of confidence in the results obtained using this new technique. 

7.1. Gas Temperature 

Knowledge of the HI gas temperature is necessary for estimating the non-thermal 
linewidths of ^^CO as well as the observed linewidths of the HI gas, the measurement of 
the HI optical depth, and column densities both for HI and molecular tracers. As such, 
accurate temperature estimates are critical. At the present time the most practical method 
for estimating the HI gas temperature in these dark cloud cores is through obtaining the 
^^CO temperature. Typically, ^^CO emission is optically thick, and thus if large scale veloc- 
ity gradients are absent, as is generally the case in dark clouds, this measurement represents 
only the surface temperature of the cloud and docs not necessarily reflect the average HI 
gas temperature in the interior of the cloud. The cores of these clouds arc shielded from 
the external UV radiation field and thus are likely to be at a lower temperature than the 
surface. While this bias introduces overestimates of the HI gas temperature by only a few 
degrees, typical core temperatures are on the order of lOK thus the effect is significant. It 
may be possible to use alternate means to estimate the temperatures of other molecular 
tracers which share similar distributions as HI within a cloud. One promising candidate 
is the inversion transitions of ammonia. Detailed comparison of this or other tracers with 
^^CO would require a better understanding of the spatial distribution of cold HI necessitating 
interferometric maps of HINSA features. 

7.2. Radiometric Noise and Baselines 

21cm spectra with high signal to noise ratios are relatively easy to obtain with short inte- 
gration times. Radiometric noise should thus not significantly affect HINSA measurements. 
Obtaining good baselines on the HI spectra is important in that any artificial slopes can 
significantly alter the derived optical depths. Since HI emission is highly variable spatially, 
ON-OFF observations are impractical. Total power ON measurents are ideal if system sta- 
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bility permits this approach to be used, as reported by (ILi &: Goldsmith II2003I ). Frequency- 
switched observations can yield sufficiently good baselines, but only if system bandpass and 
calibration are sufficiently good. A very successful combination of these techniques is avail- 
able with the recently installed Arecibo L-Band Feed Array (ALFA) in conjunction with 
the Galactic-ALFA (GALFA) spectrometer at Arecibo. Least-Squares Frequency Switching 
(LSFS) allows the RF and IF bandpasses to be separately determined so that total power 
observations produce sufficiently good baselines for HINSA observations (Krco et al. 2008, 
in preparation). 



7.3. Choice of Molecular Tracer 



Li fc Goldsmith I (120031 ) showed good correlations between the central velocities, non- 
thermal linewidths and spatial distributions of HINSA and both ^^CO and C^^O. In many 
cases, ^^CO shows a better spatial correlation, but C^^O shows a better linewidth correla- 
tion. This would not be expected if the relative abundances of HI, ^^CO, and C^^O were 
constant throughout the cloud. The variations may be a reflection of the higher critical 
density of C^^O, the density dependence of the HI/H2 conversion process, and possibly CO 
depletion onto dust grains at high densities. Thus the choice of which molecule makes the 
best template for the HINSA linewidth is complex. Based on experience, using ^^CO fits 
provides somewhat better results. As a practical matter, the ^^CO lines are much brighter 
and easier to observe than those of C^^O. As a result, we have selected ^^CO as the template 
molecule for HINSA extraction. 



7.4. Distances and Foreground Gas 



The presence of intervening HI emission sources between the observer and the HINSA 
source will ten d to wash out the absorpt ion and introduce underestimates of the HINSA 
optical depth. Ivan der Werf et al. I (119881 ) introduced an adjustment factor to compensate 
for this effect, however its precise value is poorly known. This approach was followed by 
Li fc Goldsmith I (120031 ). but these authors found it to be relatively unimportant for the 
relatively nearby clouds in the Taurus Molecular Cloud Complex. A more detailed study 
of the relation between HINSA features and distance is required when dealing with more 
distant sources. 
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7.5. Component Order 

Several velocity components are often observed in molecular clouds. It is assumed 
that there are correspondingly multiple HINSA components. Equation [1] indicates that 
when multiple overlapping HINSA components are present along the same line of sight 
with appreciable optical depths, there is significance in which component is placed in the 
background, and which is in the foreground. There is no simple way to determine the 
sequence of multiple components along the same line of sight. Since the typical optical 
depths for HINSA are on the order of 0.1, the resulting ambiguity for the HI column 
densities is below 10%, smaller than the other sources of error in our technique. 

7.6. Molecular Emission Fits and the Uniqueness Problem 

The center velocities of an individual molecular component may vary significantly within 
a cloud. It is often the case that two components may be widely separated in velocity at 
one point in the cloud while being nearly indistinguishable in another region. This poses 
a problem when trying to make measurements on individual components and is especially 
pronounced when working with HI spectra with high thermal linewidths. Further, it is 
difficult to predict how a specific implementation of the technique proposed here might react 
when working with such spectra. When working with a sufficiently small number of spectra 
it is possible to recognize and dismiss the affected spectra by eye. However, in larger regions, 
where there may be many thousands of spectra, and where fully automated fitting algorithms 
must be used, the problem becomes more pronounced. Thus we have developed a method for 
partly alleviating the uniqueness problem and allowing for easier identification of overlapping 
emission spectra. This is a problem most often avoided in many projects involving complex 
emission spectra since usually only the total integrated intensity is used. For the technique 
presented here however, the individual properties of each emission component are necessary. 

Molecular emission spectra are typically characterized by (overlapping) gaussian com- 
ponents. It is a property of gaussians that taking any gaussian function to the fourth power 
effectively reduces its linewidth by a factor of 2. Therefore, by taking all our molecular 
spectra to the fourth power prior to fitting, the individual emission components become sig- 
nificantly more apparent. However, doing this also greatly magnifies the noise in the spectra. 
Thus it is necessary to perform some filtering by removing the high-frequency Fourier com- 
ponents of each spectrum. The magnification and filtering processes can significantly alter 
the amplitude, shape, and linewidths of the final spectra. However, the final center velocities 
remain the same. Thus, the filtered spectra may be used to identify the number and center 
velocity of the emission components within a spectrum. Subsequently, the known velocities 
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can be used to constrain the fitting of tlie original data. Wliile not yielding a unique solu- 
tion, tliis metliod does allow more accurate fitting of molecular emission components using 
automated algorithms and dealing with large spectral data sets. 

8. Summary and Future Work 

The significance of accurately measuring the HI content of dark clouds and star forming 
regions has long been recognized. HINSA features have been shown as a promising method 
to achieve this goal, but the difficulty in confidently disentangling HINSA from the galactic 
background emission has limited research efforts in this field. The technique described here 
builds upon previous work to provide new opportunities. By utilizing a second derivative 
representation in which HINSA becomes the dominant feature in the spectrum, and using 
information from associated molecular tracers to constrain our fits, we are able to obtain HI 
column densities with greater confidence than previously possible. This technique enables 
us to study individual velocity components within a cloud. Several uncertainties and sources 
of error still remain, most notably the errors inherent in temperature determination through 
^^CO, and the precise relations between the properties of cold HI and molecular tracers such 
as ^^CO. As shown in f|6]the purely statistical errors derived from our constructed data were 
only a few percent in simple cases which represent the majority of observed spectra. The 
nonlinearity error may be as large as 50 percent in the extreme cases. 

While the scope of this paper is limited to a demonstration of the new technique, the 
improved confidence in the results have the potential to yield significant scientific advances 
in the field of molecular cloud and star formation studies which are deferred to a subsequent 
publication. The two astronomical sources discussed briefly here for demonstration purposes 
are part of a much larger survey of over 30 dense cloud cores. In conjunction with numerical 
modeling these data provide us with the chemical ages of the clouds and individual compo- 
nents therein, thus providing a significant constraint on theoretical models of star formation. 
Large maps have been collected with the Arecibo L-Band Feed Array (ALFA) of the Taurus 
and Perseus star forming regions. These regions show plentiful HINSA features whose anal- 
ysis may shed light on the dynamics of such complexes including the processes which may 
have triggered their formation. 
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A. Description of the Procedure 

For definitiveness, we here describe the process of obtaining molecular parameters to be 
used as a template as well as the actual HINSA extraction in a step-by-step fashion. We 
model this description based on our specific observations and available data although the 
technique can be easily adapted to other circumstances. We obtained ^^CO and ^^CO data 
with the 14 m Five College Radio Observatory having a 45" to 50" FWHM beam size and 
the HI data were obtained using the 100 m Green Bank Telescope having a 9' FWHM beam 
size. 



1. Reduction and Regridding: It is assumed that all data have been properly reduced, 
calibrated, and regridded. 

2. Convolution: Initially the ^^CO and ^^CO spectra are convolved to a 2' FWHM beam 
size in order to improve the signal to noise ratio. A gaussian convolving beam shape 
is used. 

3. CO Fitting (first pass): While assuming a particular gas temperature (lOK), emission 
functions are fitted to all ^^CO spectra to obtain their center velocities. These velocities 
are used to determine the excitation temperature of each ^'^CO velocity component, 
using the accompanying ^^CO spectrum. The temperatures are derived under the 
assumpti on that the CO e missi on is in LTE and optically thick at line center. Then, 



for i^CO, IStahler fc PaUa I (120041 ) 



/(T,,„) = /(T,.) = ^ + fin,) , (Al) 

where the excitation temperature (Tex) is taken to be equal to the kinetic temperature 
(Tkin), is the observed brightness temperature at line center. To is the equivalent 
temperature of the transition (5.5K for the J=l - transition of ^^CO), T;,^ is the 
blackbody radiation temperature of the background field (2.7K), and 

fiT) = ^-^r^ . (A2) 

■'^ ' exp{To/T)-l ^ ' 

4. CO fitting (second pass): Linewidth, and optical depth are derived for each ^^CO 
emission component using the previously obtained temperatures and center velocities. 
Subsequently the non-thermal linewidths are obtained using: 

o^L = (^It + , (A3) 
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where aobs is the observed total hnewidth, and ant and ath are respectively the non- 
thermal and thermal components of the linewidth. 

5. Molecular Column Densities: Column densities for ^'^CO for each component along 



each l ine of sight are calculated according to the technique described in lStahler &: Palla 



( 120041 ) ■ H2 column densities are estimated using a H2/ ^^CO abundance ratio of 7.5 x 10^ 



for these galactic clouds. 

6. Estimating Atomic Parameters: The center velocity, non-thermal linewidth, and tem- 
perature for the HI gas are determined for each velocity component along each line of 
sight by convolving the previously obtained ^^CO parameters to the GET 9' beam as 
weighted by the ^^CO column densities according to: 



CO 
i,comp 



-y /-^t •) \ i-,J / - - i.,i^uiiLiJ - i,comp { AA") 



^CO 

* J y ' i,comp 



where Xj^comp represents the temperature, non-thermal linewidth, or center velocity of 
a particular velocity component at sky position j, i represents all nearby sky positions 
to be convolved, f{Oij) represents the beam response function for angular distance 9 
between points i and j, and A^i^comp represents the corresponding ^^CO column density 
for a particular component at position i. For comparison with HI column densities, 
the H2 column densities are similarly convolved without the mass weighting factor. 

7. Having estimates of the temperature, central velocity, and linewidth for each HINSA 
component we can, for any value of the optical depth, remove the HINSA absorption 
to obtain the background emission using equation [H Since their true ordering along 
the line of sight is not known, as a matter of definition velocity components with the 
highest positive velocities are taken to be the most distant. 



8. Finally we search for those values of Tq and an which minimize the value of equation 
[H]to arrive at our solutions. 







HI Mass(M0) 


H2 Mass(MQ) 


Mean HI/H2 


Peak H2 Dens, (cm 3) 


Mean H2 Dens. (cm 3) 


L134 A 


2.35 - 3.40 


4.0 • 10-2 


71 


1.1 • 10-3 


2490 


910 


L134 B 


1.80-2.75 


9.5 • 10-3 


41 


4.7-10-^ 


2390 


660 


L134 C 


0.50 - 2.75 


4.7- 10-3 


8.4 


1.1 • 10-3 


1530 


310 


L134 D 


-0.7-1.10 


5.1 • 10-3 


12 


8.5 • 10-^ 


1230 


480 


L1757 A 


4.70 - 6.00 


1.6 ■ 10-2 


31 


1.0 • 10-3 


840 


310 


L1757 B 


2.50 - 3.35 


2.7-10-^ 


45 


1.1 ■ 10-2 


540 


270 


L1757 C 


1.50-2.20 


3.0 ■ 10-1 


106 


5.7-10-3 


1560 


305 


L1757 D 


0.80 - 1.50 


1.7-10-1 


95 


3.5 • 10-3 


880 


290 



Table 1: Derived Component Properties for L134 and L1757. The component radii are estimated by the extent of 
the i3C0 emission and are used to estimate cloud depths and thus volume densities. In L1757, the fact that some 
components have roughly the same H2 density but very different HI abundances suggests that the different components 
may have formed at different times. 





Peak HINSA optical depth 


Peak ^^CO optical depth 


Peak i^CO Tex (K) 


Mean continuum at 21cm (K) 


L134 A 


0.53 


0.84 


19 


3.1 


L134 B 


0.20 


0.67 


17 




L134 C 


0.20 


0.24 


17 




L134 D 


0.19 


0.22 


14 




L1757 A 


0.20 


0.31 


16 


3.2 


L1757 B 


0.35 


0.29 


17 




L1757 C 


0.36 


0.74 


19 




L1757 D 


0.30 


0.35 


15 





Table 2: Derived Component Properties for L134 and L1757 continued. The CO temperature, derived directly from 
the ^^CO spectra, is used as an estimate of the HI temperature represented by in Equation [TJ 
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Input Optical Depth 0,2 
Residual Minimum at 0.145 




0.0 0.1 0.2 0.3 0.4 

Trial Peak Optical Depth 



Fig. 12. — A similar residuals plot as in Figure [H] as applied to the solution in Figure [TTl 
The behavior of the residuals is the same here even though the error in the solution is clearly 
evident (r = 0.15 compared to 0.20). This demonstrates the difficulty in identifying incorrect 
fits in cases such as in Figure [TTl 



